AFRL-RX-WP-TP-2008-4354 


MODELING  OF  TEXTURE  EVOLUTION  DURING  HOT 
FORGING  OF  ALPHA/BETA  TITANIUM  ALLOYS 
(PREPRINT) 


M.G.  Glavicic,  R.L.  Goetz,  D.R.  Barker,  G.  Shen,  D.  Furrer,  A.  Woodfield,  and  S.L.  Semiatin 


Metals  Branch 

Metals,  Ceramics,  and  NDE  Division 


JUNE  2007 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
MATERIALS  AND  MANUFACTURING  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7750 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 

Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY) 

June  2007 

2.  REPORT  TYPE 

Journal  Article  Preprint 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

MODELING  OF  TEXTURE  EVOLUTION  DURING  HOT  FORGING  OF 
ALPHA/BETA  TITANIUM  ALLOYS  (PREPRINT) 

5a.  CONTRACT  NUMBER 

In-house 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

62102F 

6.  AUTHOR(S) 

M.G.  Glavicic,  R.L.  Goetz,  and  D.R.  Barker  (UES,  Inc.) 

G.  Shen  and  D.  Furrer  (Rolls-Royce  Corporation) 

A.  Woodfield  (General  Electric  Aviation) 

S.L.  Semiatin  (AFRL/RXLMP) 

5d.  PROJECT  NUMBER 

4347 

5e.  TASK  NUMBER 

RG 

5f.  WORK  UNIT  NUMBER 

M02R2000 

7.  PERFORMING  ORGANIZATION  NAME) 

UES,  Inc. 

S)  AND  ADDRESS(ES) 

Metals  Branch  (AFRL/RXLMP) 

Metals,  Ceramics,  and  NDE  Division 

Materials  and  Manufacturing  Directorate 
Wright-Patterson  Air  Force  Base,  OH  45433-7750 

Air  Force  Materiel  Command,  United  States  Air  Force 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-RX-WP-TP-2008-4354 

Rolls-Royce  Corporation 

General  Electric  Aviation 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

Materials  and  Manufacturing  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7750 

Air  Force  Materiel  Command 

United  States  Air  Force 

10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/RXLMP 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-RX-WP-TP-2008-4354 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

Journal  article  submitted  to  Metallurgical  and  Materials  Transactions  A.  PAO  Case  Number:  AFRL/WS  07-1466; 
Clearance  Date:  20  Jun  2007.  The  U.S.  Government  is  joint  author  of  this  work  and  has  the  right  to  use,  modify, 
reproduce,  release,  perform,  display,  or  disclose  the  work. 

14.  ABSTRACT 

Texture  development  during  the  subtransus  hot  forging  of  alpha/beta  titanium  alloys  with  an  equiaxed-alpha 
microstructure  was  modeled  using  the  Los  Alamos  Polycrystalline  Plasticity  (LApp)  code  and  <1 1 1>  pencil-glide 
polycrystalline  plasticity  code  coupled  with  the  finite-element-method  (FEM)  program  DEFORM.  The  methodology 
treated  the  partitioning  of  the  imposed  strain  between  the  alpha  and  beta  phases,  and  thus  enabled  the  prediction  of  the 
distinct  deformation  textures  developed  in  the  primary  alpha  and  beta  matrix  during  hot  working.  Two  variant  selection 
rules  in  conjunction  with  the  beta  deformation  texture  were  also  examined  to  establish  a  method  for  predicting  the 
transformation  texture  of  secondary  alpha  developed  as  a  result  of  beta  decomposition  during  cool-down  following 
forging  or  heat  treatment.  The  approach  was  validated  via  an  industrial-scale  trail  comprising  hot  pancake  forging  of  Ti- 
6A1-4V. 


15.  SUBJECT  TERMS 

titanium,  texture,  modeling,  strain  partitioning,  variant  selection 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT: 

SAR 

18.  NUMBER 

OF  PAGES 

30 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

Sheldon  L.  Semiatin 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


1 


Modeling  of  Texture  Evolution  during  Hot  Forging  of  Alpha/Beta 

Titanium  Alloys 

M.G.  Glavicic*,  R.L.  Goetz*,  D.R.  Barker*,  G.  Shen**, 

D.  Furrer**,  A.  Woodfleld***,  and  S.L.  Semiatin 
Air  Force  Research  Laboratory,  Materials  and  Manufacturing  Directorate, 
AFRL/MLLM,  Wright-Patterson  AFB,  OH  45433-7817 
*UES,  Inc.,  4401  Dayton-Xenia  Road,  Dayton,  OH  45432 
**Rolls-Royce  Corporation,  P.O.  Box  420.  Indianapolis,  IN  46206-0420 
***General  Electric  Aviation,  Cincinnati,  OH  45215 
Abstract 

Texture  development  during  the  subtransus  hot  forging  of  alpha/beta  titanium 
alloys  with  an  equiaxed-alpha  microstructure  was  modeled  using  the  Los  Alamos 
Poly  crystalline  Plasticity  (LApp)  code  and  a  <1 1 1>  pencil-glide  poly  crystalline  plasticity 
code  coupled  with  the  finite-element-method  (FEM)  program  DEFORM™.  The 
methodology  treated  the  partitioning  of  the  imposed  strain  between  the  alpha  and  beta 
phases,  and  thus  enabled  the  prediction  of  the  distinct  deformation  textures  developed  in 
the  primary  alpha  and  beta  matrix  during  hot  working.  Two  variant  selection  rules  in 
conjunction  with  the  beta  deformation  texture  were  also  examined  to  establish  a  method 
for  predicting  the  transformation  texture  of  secondary  alpha  developed  as  a  result  of  beta 
decomposition  during  cool-down  following  forging  or  heat  treatment.  The  approach  was 
validated  via  an  industrial-scale  trial  comprising  hot  pancake  forging  of  Ti-6A1-4V. 
Keywords:  Titanium,  Texture,  Modeling,  Strain  Partitioning,  Variant  Selection 


1.  Introduction 

The  modeling  of  texture  evolution  during  the  hot  working  of  two-phase, 
alpha/beta  titanium  alloys  with  a  microstructure  of  equiaxed  (primary)  alpha  in  a  matrix 
of  beta  is  complicated  by  several  factors.  First,  the  flow  stress  behaviors  of  the 
hexagonal-close-packed  (hep)  alpha  phase  and  the  body-centered-cubic  (bee)  beta  phase 
differ  and  exhibit  different  dependences  on  temperature.  Such  differences  result  in  an 
unequal  partitioning  of  the  imposed  strain  that  drives  the  formation  of  deformation 
texture.  Second,  the  beta  phase  decomposes  to  form  secondary  (platelet)  alpha  during 
cooling  from  the  hot  working  (or  final  heat  treatment)  temperature.  At  moderate  cooling 
rates,  this  phase  transformation  follows  a  Burgers-type  orientation  relationship, 

{HO}^  //(0001)a  (1) 

<111>3  //[2ll0]„ 

in  which  there  are  twelve  distinct  possible  variants  that  can  form  from  a  single  orientation 
of  a  prior  beta-phase  grain;  i.e.,  the  secondary-alpha  orientation  is  related  to  one  of  six 
{ 1 10}p  planes,  each  of  which  contains  two  <11  l>p  directions.  In  addition,  the 
transformation  texture  thus  formed  can  vary  widely  depending  upon  the  initial  texture  of 
the  beta-phase  and  any  predisposition  for  one  or  several  of  the  twelve  variants  to  form 
preferentially  over  the  others. 

The  objective  of  the  present  work  was  to  develop  a  modeling  framework  to 
describe  the  evolution  of  deformation  and  transformation  texture  during  the 
thermomechanical  processing  of  alpha/beta  titanium  alloys  at  temperatures  in  the  two- 
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phase  field.  For  this  purpose,  numerical  techniques  to  describe  strain  partitioning  and  the 
crystal  rotations  associated  with  metal  flow  and  crystal  plasticity  (i.e.,  slip)  were 
integrated  in  a  user-friendly  manner  and  validated  using  an  industrial-scale  forging  of  Ti- 
6A1-4V. 

2.  Background 

A  brief  review  of  pertinent  previous  efforts  in  the  area  of  texture  modeling  is 
presented  below  to  provide  a  background  for  the  methodology  used  in  the  present  work. 

2.1  Texture  Modeling 

Several  types  of  models  can  be  used  to  predict  the  evolution  of  crystallographic 
texture  during  the  large  strain  deformation  of  polycrystalline  materials.  These  include  the 
Sachs  isostress  approach  [1],  the  Taylor  isostrain  model  [2],  and  various  crystal  plasticity 
fmite-element-method  (CPFEM)  codes.  The  Taylor  model,  which  enforces  strain 
compatibility  between  grains,  generally  provides  the  best  balance  between  simplicity  and 
accuracy.  It  forms  the  basis  for  the  widely-used  Los  Alamos  Polycrystalline  Plasticity 
(LApp)  code  [3],  which  can  be  used  to  establish  the  evolution  of  crystallographic  texture 
in  single-phase  metals  during  rolling,  forging,  extrusion,  etc.  However,  the  Taylor 
approach  does  not  allow  for  strain  variations  from  grain-to-grain,  let  alone  within 
individual  grains.  When  such  non-uniformity  is  important,  viscoplastic  CPFEM  models 
[4-6]  are  applied,  but  require  substantial  computer  resources  and  are  significantly  more 
complex  to  apply  to  the  large-strain  deformation  involved  in  an  actual  forging  process. 
Recently,  the  program  ALAMEL  [7]  was  developed  in  an  attempt  to  reduce  computation 
requirements;  it  was  applied  to  model  texture  evolution  during  the  rolling  of  interstitial- 
free  steel.  In  other  work,  Lebensohn  and  Canova  [8]  developed  a  viscoplastic  self- 
consistent  model  to  describe  texture  evolution  in  alpha/beta  titanium  alloys  with  a 
colony-alpha  microstructure.  Subsequently,  this  approach  was  validated  for  the  simple 
compression  and  rolling  of  the  two-phase  Fe-Cu  system  by  Commentz,  et  al.  [9]. 
Although  each  of  these  models  shows  promise,  neither  is  yet  suitable  for  quantifying 
texture  evolution  during  the  hot  forging  of  complex  two-phase  alloys  in  an  industrial 
environment. 

2.2  Modeling  of  Metal  Flow  during  Forging 

Rigid,  viscoplastic  FEM  codes,  such  as  DEFORM™  [10],  are  now  used 
commonly  in  the  forging  industry.  Such  software  has  been  applied  to  a  variety  of  non¬ 
steady  and  steady-state  processes  to  predict  the  effect  of  material  properties,  die  design, 
and  process  variables  on  metal  flow  and  the  occurrence  of  defects.  Typically,  the 
workpiece  material  is  assumed  to  be  an  isotropic  continuum  (i.e.,  with  a  random  texture) 
whose  deformation  behavior  is  described  by  the  von  Mises  criterion  and  associated  flow 
rule.  As  a  consequence,  the  development  of  crystallographic  texture  is  not  taken  into 
account,  let  alone  the  effect  of  such  texture  evolution  on  the  metal  flow  itself. 

2.3.  Coupled  Metal  Flow-Texture  Models 

To  simplify  the  problem  associated  with  describing  metal  flow  and  texture 
evolution  during  large-strain  deformation,  several  research  efforts  have  sought  to  couple 
continuum  FEM  and  crystal-plasticity  codes.  Perocheau,  et  al.  [11]  coupled  an  FEM 
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code  with  a  polycrystalline  plasticity  model  based  upon  the  relaxed-constraints  (RC) 
Taylor  hypothesis  to  model  the  extrusion  and  rolling  of  fee  metals.  Schoenfeld,  et  al. 
[12],  modeled  the  through-thickness  rolling  textures  of  aluminum  by  including  advanced 
friction  and  internal-variable  constitutive  models  into  a  2D  FEM  model  developed  by 
Lee,  et  al.  [13].  More  recently,  the  ALAMEL  model  [7],  a  statistical-type  method  with 
Taylor  type  homogenization,  was  used  to  model  the  texture  developed  in  rolled 
interstitial-free  steel. 


3.  Modeling  Approach 

In  a  similar  fashion  to  the  previous  efforts  of  Van  Houtte,  et  al.,  Perocheau,  et  al., 
and  Schoenfeld,  et  al.  [7,  11,  12]  on  single -phase  alloys,  the  present  work  to  analyze 
texture  evolution  during  the  hot  forging  of  two-phase  alpha/beta  titanium  alloys,  used  a 
two  dimensional  FEM  program  to  analyze  crystal  rotations  due  to  metal  flow  and  two 
Taylor  isostrain  models  (Los  Alamos  Poly  crystalline  Plasticity  (LApp)  code  [3]  for  the 
alpha  phase  and  a  <1 1 1>  pencil-glide  poly  crystalline  plasticity  code  for  the  beta  phase)  to 
quantify  rotations  due  to  crystal  plasticity.  The  distinguishing  features  of  the  present 
work  related  to  the  treatment  of  strain  partitioning  between  the  phases  and  its  effect  on 
deformation  texture  evolution,  the  description  of  transformation  texture  and  variant 
selection  rules,  and  the  integration  of  the  software  packages. 


3.1  Self-Consistent  Model  of  Strain  Partitioning 

The  flow  stresses  of  the  individual  phases  in  alpha/beta  titanium  alloys  such  as  Ti- 
6A1-4V  are  very  different  at  the  temperatures  typically  used  during  hot  forging,  i.e., 
800°C  to  975°C.  The  alpha  phase  is  approximately  three  times  as  strong  as  the  beta 
phase  [14],  In  addition,  the  volume  fraction  of  the  phases  varies  greatly.  For  example, 
the  volume  fraction  of  alpha  varies  from  ~80  pet.  at  800°C  to  ~20  pet.  at  975°C.  Hence, 
the  strain  accommodated  by  the  alpha  and  beta  phases  is  not  equal  and  may  be  expected 
to  vary  with  temperature.  Thus,  an  analytical  method  to  quantify  the  strain 
accommodated  by  each  phase  is  required  to  model  the  texture  change  associated  with 
crystal  plasticity.  In  the  present  work,  a  simple  self-consistent  approach  [14]  was 
adopted.  This  method  is  based  on  the  technique  developed  by  Hill  [15]  and  later 
extended  by  Suquet  [16]  for  linearly-elastic  solids.  Hill’s  analysis  assumed  that  both 
phases  are  linearly  viscoplastic;  i.e.,  they  each  have  a  constitutive  relation  of  the  form: 


<J,  =  k; 


(2) 


in  which  a  and  e  denote  the  flow  stress  and  strain  rate,  respectively,  kL  is  the  “viscosity” 
coefficient,  and  the  subscripts  (i  =  1,  2)  refer  to  phases  1  or  2.  The  viscosity  of  the 


aggregate  k^c  (which  relates  the  aggregate  flow  stress  and  strain  rate)  is  given  by  the 


following  expression: 

k^/k;=(l/6)j3-2j0+5(l-/)(j0-l)  +  ^[3-2p  +  5(l-/)(p-l)]2  +  24pj,(3) 


in  which  p  =  kj/kJ- ,  and  f  denotes  the  volume  fraction  of  phase  1. 

Suquet  [16]  extended  the  above  analysis  to  the  case  in  which  both  phases  are 
power-law  viscoplastic,  viz., 


cr.  =k.£" 


(4) 
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and  the  strain-rate  sensitivity  exponents  for  both  phases  are  equal;  i.e.,  mi  =  m2  =  m.  In 
this  case,  the  viscosity-like  parameter  of  the  aggregate  (k)  is  a  function  of  the  values  of 
viscosity-like  parameters  for  the  two  phases,  ki  and  k2,  as  well  as  m,  p,  and  ksLc/k^ : 


k  /  kx 


min 

p  >  0 


(m+l)/2 


(m+ 1)/ 

f  +  (l-f)p 


(5) 


In  Equation  (5),  the  expression  on  the  right-hand  side  before  the  argument  in  braces 
denotes  the  minimum  value  for  p  >  0.  In  practice,  the  value  of  ksLc/k,L  for  the 
corresponding  linearly-viscous  case  is  unknown.  Thus,  the  evaluation  of  k/ki  involves  a 
procedure  in  which  trial  values  of  p  are  guessed,  k^./kf',  is  determined  from  Equation 

(3),  and  p  and  k^/kf-  are  inserted  into  Equation  (5).  The  value  of  p  that  yields  the 

minimum  value  of  k/ki  is  the  appropriate  one. 

The  average  strain  rates  in  the  two  phases  are  readily  calculated  from  the  values 
of  k,  ki,  k2,  and  the  volume  fraction  f  of  phase  “1”.  Following  Hill  [17],  the  aggregate 
flow  stress  gov  and  strain  rate  sm  are  volume  averages  of  the  corresponding  flow  stresses 
and  strain  rates  in  the  individual  phases: 

CTov  =  =  fV” + (l  -f )  M2”  >  (6) 

=  fri+(i-f)4  »  (?) 

Solving  Equation  (7)  for  s2  as  a  function  of  sov  and  ex , 

e2leow=[\-f(sJeow)]l(\-f).  (8) 

Inserting  this  relation  into  Equation  (6),  an  expression  for  sx  /  eov  is  obtained: 


Equation  (9)  cannot  be  solved  analytically,  but  is  readily  evaluated  using  numerical 
techniques. 

From  Equations  (8)  and  (9),  the  strain  rates  and  hence  strain  increments  for  each 
phase  are  determined  as  a  function  of  their  relative  flow  stresses  and  volume  fractions. 
The  strain  increments  so  determined  were  used  in  crystal-plasticity-model  calculations  to 
estimate  rotations  due  to  slip. 


3.2  Beta-Phase  Decomposition 

At  relatively  slow  cooling  rates  following  hot  forging  (or  heat  treatment)  of 
titanium  components  (i.e.  ,  <  20°C/min),  primary  alpha  particles  grow  epitaxially, 
consuming  the  majority  of  the  beta  matrix  [8],  At  faster  cooling  rates  (between  ~20  and 
200°C/min),  typical  of  fan  cooling  or  water  quenching  of  finite-section-size  workpieces, 
the  beta  matrix  decomposes  at  a  critical  supersaturation  to  form  colonies  of  secondary- 
alpha  platelets  [18]  each  of  which  possess  a  Burgers-orientation  relation  with  its  parent 
beta  matrix  (Equation  (1)). 

The  texture  associated  with  the  secondary-alpha  phase  can  vary  greatly  depending 
on  which  of  the  twelve  possible  variants  is  selected  during  beta-phase  decomposition. 
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The  manner  in  which  a  beta  phase  orientation  is  transformed  is  governed  by  the 
relationship 

DSPB  =  A,  (10) 

in  which  D  is  a  rotation  matrix,  defined  by  the  Euler  angles  (pi  =  135,  ®  =  90,  cp2  =  325 
[19],  is  a  subset  of  the  cubic  symmetry  operators  (Table  1),  B  is  the  beta  phase 

orientation  matrix,  and  A  is  the  alpha  phase  orientation  matrix. 

In  the  present  work,  the  secondary  alpha  of  interest  was  formed  during  rapid 
cooling  following  solution  treatment  of  a  forging  high  in  the  alpha-beta  phase  field.  Two 
different  variant-selection  possibilities  were  examined.  These  rules  were  based  on  the 
previous  experimental  observations  for  commercial-purity  titanium  [20-22]  and  Ti-6A1- 
4V  [23-25], 

In  the  first  case,  all  variants  were  assumed  to  be  equally  likely.  In  this  instance,  a 
specific  rotation  matrix  Sk  from  Table  1  was  randomly  selected  from  the  twelve 

possibilities  for  the  transformation  of  each  individual  beta  phase  orientation. 

In  the  second  case,  <1 1 1>  pencil-glide  in  the  beta  phase  during  cooling  from  the 
final  solution  heat  treatment  was  assumed  to  control  variant  selection.  Here,  the 
methodology  developed  by  Morris  and  Semiatin  [26]  was  used  to  determine  the  <11 1> 
slip  directions  and  slip  planes  activated  to  accommodate  the  strain  increments  developed 
during  cooling,  as  predicted  by  FEM  analysis.  From  the  two  closest  {110}  planes  that 
could  have  equivalently  accommodated  this  same  shear  strain  via  <1 1 1>  pencil-glide  [see 
Figure  4  in  reference  34],  the  { 1 10}<1 1 1>  slip  system  on  which  the  highest  resolved 
shear  stress  had  been  developed  was  assumed  to  comprise  the  plane/direction  that 
determined  the  orientation  of  the  selected  alpha  variant.  Because  the  variants  of  the  alpha 
phase  are  assumed  to  form  per  the  Burgers  relationship  (Equation  (1)),  the  resolved  shear 
stress  x  on  a  given  {1 10}  <  1 1 1  >  slip  system  in  the  beta  phase  for  a  given  imposed  stress 

state  will  be  the  same  as  the  resolved  shear  stress  on  the  (0001)[2  110]  slip  system  in  the 
alpha  phase  after  the  phase  transformation.  The  equivalent  slip  planes/directions  of  the 
beta  and  alpha  phase  are  then  related  by 

t=  mf  :BcrBT  =  mb“sal<a>  ;  A  cr  A1 ,  (11) 

in  which  m(jf  and  mbaasal  <a>  are  the  tensor  quantities  m\.  )  =  bfi  +  bjni  ,  in  which  the 

slip  direction  b ,  and  slip  plane  normal  h  are  defined  by  {1 10}  <  1 1 1  >  slip  systems  of  the 
beta  phase  and  the  (0001)  [2  1  1 0]  slip  system  of  the  alpha  phase.  Inserting  Equation  (10) 
into  Equation  (11) 

mf  :  B  <y  BT  =  mbaasal<a>  :  D  SPB  cr  (d  Sp  bJ  ,  (12) 

permits  the  specific  alpha  phase  variant  D  B  that  has  the  same  resolved  shear  stress 

as  in  the  beta  phase  to  be  selected  from  the  12  possible  variants  and  used  to  transform  the 
beta  phase  orientation. 

3.3  Calculation  Procedures  and  Integration  of  Software  Packages 

Current  crystal-plasticity  FEM  codes  are  too  demanding  from  a  computational 
standpoint  to  model  texture  evolution  at  every  material  point  during  the  forging  of  a 
single-phase  alloy,  let  alone  for  a  complex  forging  of  a  two-phase  material.  Thus,  the 
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present  work  consisted  of  a  decoupled  approach  in  which  crystal  rotations  due  to  metal 
flow  and  crystallographic  slip  were  treated  separately,  thereby  providing  a  compromise 
between  scientific  rigor  and  user-friendliness  in  an  industrial  environment.  To  this  end, 
the  FEM  metal- flow  program  DEFORM™  and  the  two  crystal-plasticity  codes  (LApp  [3] 
and  the  <11 1>  pencil-glide  model  of  Morris  and  Semiatin  [26])  were  selected.  In 
addition  to  these  software  packages,  other  specially-developed  codes  were  incorporated 
into  a  series  of  user-friendly  GUIs  that  could  be  easily  operated. 

Model  calculations  consisted  of  running  DEFORM™  in  the  usual  manner.  After 
the  DEFORM™  simulation  was  run,  the  strain  increments  and  metal-flow-related  rotation 
for  each  time  step  and  FEM  node  (integration  point)  was  output  in  the  form  of  a  set  of 
ASCII  files.  The  strain  increments  were  partitioned  per  the  approach  described  in 
Section  3.1  and  used  as  input  to  separate  alpha-phase  and  beta-phase  calculations  of  the 
slip-related  rotations.  For  the  alpha-phase,  LApp  [3]  was  used  to  model  the  texture 
evolution,  whereas  for  the  beta  phase,  the  Morris  and  Semiatin  <11 1>  pencil-glide 
approach  [26]  was  used.  The  crystal  plasticity  input  also  included  experimental  pole- 
figure  data  (from  x-ray  or  neutron  diffraction  or  EBSD  measurements)  obtained  using  the 
Los  Alamos  Preferred  Orientation  Package  (popLA)  [3]  and  other  internally  developed 
software.  After  each  strain  increment,  the  texture  change  (equal  to  the  sum  of  the  metal- 
flow  and  LApp/<l  1 1> pencil-glide  calculated,  slip-related  rotations)  was  imposed  at  each 
node  for  each  phase.  The  resulting  texture  was  then  used  as  the  initial  texture  for  the  next 
increment  of  deformation. 

Subsequent  to  the  crystal  plasticity  calculations,  an  additional  FEM  simulation 
was  performed  to  model  the  quenching  procedure  that  followed  the  final  heat  treatment 
of  a  pancake  forging  used  as  the  validation  example.  The  predicted  beta-phase  texture 
was  then  converted  to  a  secondary-alpha  texture  using  each  of  the  two  possible  variant- 
selection  rules  using  the  stresses  and  strains  predicted  to  have  been  developed  at  the 
median  temperature  (-810  °C)  at  which  secondary  alpha  would  have  been  expected  to 
nucleate  during  quenching  [27],  The  predicted  primary-alpha  (deformation)  texture  and 
secondary-alpha  (transformation)  texture  were  then  weighted  per  the  relative  volume 
fractions  of  the  phases  and  added  to  obtain  the  overall  alpha-phase  texture.  The  volume 
fraction  of  primary  alpha  is  function  of  the  forging  temperature  and  the  cooling  rate;  it 
can  be  estimated  using  a  simple  diffusion  analysis  [18,  27].  The  volume  fraction  of 
secondary  alpha  was  then  taken  to  be  the  difference  between  the  total  alpha  content 
(typically  of  the  order  of  0.95  for  alloys  such  as  Ti-6A1-4V)  and  the  primary-alpha 
volume  fraction.  Alternatively,  EBSD  techniques,  such  as  that  described  by  Germain,  et 
al.  [28],  can  be  used  to  separate  the  primary-  and  secondary-alpha  constituents. 

4.  Materials  and  Experimental  Procedures 

The  texture  developed  during  the  final  heat  treatment  of  a  Ti-6A1-4V  pancake 
forging  was  used  to  validate  the  modeling  approach  described  in  Section  3. 

4.1  Material 

The  validation  material  comprised  a  280-mm  diameter  Ti-6A1-4V  billet  whose 
composition  (in  weight  percent)  was  6.62  aluminum,  4.10  vanadium,  0.19  oxygen,  0.007 
nitrogen,  0.02  carbon,  0.004  hydrogen,  balance  titanium.  The  initial  microstructure  of  the 
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billet  material  at  the  forging  preheat  temperature  (955°C)  consisted  of  equiaxed,  primary- 
alpha  particles  in  a  matrix  of  the  beta  phase  (Fig.  1). 

The  variation  of  the  volume  fraction  of  beta  with  temperature  (the  so-called  beta- 
approach  curve),  which  was  needed  for  both  the  strain-partitioning  calculations  and  the 
determination  of  the  overall  alpha-phase  texture,  was  determined  using  PANDAT™  [29] 
and  its  associated  titanium  database.  The  beta-approach  curve  between  700°C  and  the 
beta-transus  temperature  (~995°C)  was  fit  to  a  seventh-order  polynomial  function. 

4.2  Experimental  Procedures 

The  validation  trial  comprised  the  hot-die  forging  of  a  cylindrical  billet  whose 
initial  dimensions  were  280-mm  diameter  x  510-mm  height.  The  preform  was  coated 
with  glass,  preheated  at  955°C,  and  forged  to  a  70-pct.  height  reduction  between  flat  tool- 
steel  dies  preheated  to  approximately  600°C  to  minimize  workpiece  chilling.  Forging 
was  conducted  in  a  hydraulic  press  at  a  constant  crosshead  speed  to  produce  a  nominal 
strain  rate  of  0.7  s'1.  With  the  prescribed  die  and  billet  temperatures  and  deformation 
rate,  it  was  estimated  that  the  preform  temperature  remained  between  870  and  955°C  and 
that  deformation  heating  at  the  center  of  the  preform  did  not  play  a  significant  role. 
Immediately  following  forging,  the  pancake  was  water  quenched.  Subsequently,  the 
forging  was  given  a  final  heat  treatment  at  a  temperature  30°C  below  the  beta  transus  for 
one  hour  and  water  quenched. 

The  forging  and  quenching  operations  were  modeled  using  DEFORM™.  The 
point-tracking  feature  within  the  software  was  used  to  identify  the  local  strain  increments, 
metal-flow  rotation,  and  cooling  transients  (during  water  quenching)  at  six  specific 
locations  within  the  pancake  (Figures  2  and  3). 

Specimens  extracted  from  the  pancake  forging  at  the  six  tracking-point  locations 
were  metallographically  prepared  for  microstructural  characterization  and  texture 
measurement.  The  microstructures  were  characterized  using  a  Leica  Stereoscan  360  FE 
scanning-electron  microscope  (SEM).  The  textures  in  the  final  forging,  as  well  as  at  the 
corresponding  locations  in  the  starting  billet,  were  quantified  using  x-ray  diffraction  to 
establish  (1010),  (0002),  (101  1),  (1012),  and  (1 120)  alpha-phase  pole  figures  and 
(110)  and  (200)  beta-phase  pole  figures.  The  pole  figures  were  measured  in  the  Schultz- 
reflection  mode  using  Cu  Ka  radiation.  Analysis  of  the  measured  data  was  performed 
using  popLA.  The  input  to  crystal-plasticity  codes  for  the  deformation-texture 
simulations  consisted  of  the  initial  textures  and  slip  modes.  To  model  the  evolution  of 
texture  in  the  primary-alpha  phase,  LApp  [3]  was  used  and  deformation  was  assumed  to 
be  possible  via  basal  <a>,  prism  <a>,  and  pyramidal  <c+a>  slip  with  a  ratio  of  the 
corresponding  critical-resolved  shear  stresses  of  1,  0.7,  and  3,  respectively.  For  the  beta 
phase,  deformation  via  <11 1>  pencil-glide  was  modeled  using  the  methodology  of 
Morris  and  Semiatin  [26]. 

In  order  to  gage  the  accuracy  of  the  modeling  capability  for  the  alpha-phase 
deformation  texture  alone,  a  specimen  from  the  RT  location  in  the  forging  was  taken  to 
the  final  heat  treatment  temperature  and  then  furnace  cooled  to  produce  a  structure  that 
was  essentially  fully  equiaxed  alpha.  Pole  figure  measurements  for  this  sample  were 
made  in  the  same  manner  as  the  other  samples. 
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5.  Results  and  Discussion 


5.1  Microstructure  Observations 

Following  the  final  heat  treatment,  a  variety  of  bimodal  microstructures  was 
developed  in  the  forging  (Fig.  4).  The  different  microstructures  can  be  ascribed  to  the 
variation  in  local  cooling  rate  (Fig.  3),  the  supersaturations  that  were  developed,  and  the 
temperature  at  which  the  secondary  alpha  was  formed  [18].  At  the  bore-midplane 
location  (BC,  Fig.  4a),  for  example,  a  bimodal  microstructure  of  equiaxed  alpha  in  a 
transformed  matrix  of  coarse,  colony  (secondary)  alpha  was  formed.  The  microstructure 
was  similar  at  the  top  of  the  bore  (location  BT,  Fig.  4c),  despite  the  somewhat  faster 
cooling  rate  at  BT  compared  to  that  at  BC  (Fig.  3).  At  the  rim-midplane  location  (RC, 
Fig.  4b),  the  secondary  alpha  was  somewhat  finer  even  though  the  cooling  rate  here  (Fig. 
3)  was  estimated  to  be  between  that  of  the  bore-midplane  and  bore-top  locations.  A 
much  finer  microstructure  was  formed  at  the  rim-top  region  (RT,  Fig.  4d);  higher 
magnification  micrographs  revealed  that  the  secondary  alpha  here  had  a  basketweave 
morphology  whose  formation  could  be  explained  by  the  very  fast  cooling  rate,  deduced 
from  the  FEM  simulations  (Fig.  3).  In  general,  the  variation  of  the  volume  fraction  of 
primary  alpha  with  location  (Table  2)  corresponded  to  the  cooling  rate;  i.e.,  slower 
cooling  rates  led  to  a  greater  tendency  for  epitaxial  growth  of  primary  alpha  and  hence  a 
greater  final  fraction  of  this  constituent.  On  the  other  hand,  a  more  detailed  analysis 
appears  to  be  needed  to  explain  the  details  of  the  formation  and  thickening  of  secondary- 
alpha  plates. 

5.2.  Texture  Measurements  and  Model  Predictions 

The  texture  developed  in  the  pancake  forging  also  varied  from  location  to 
location.  At  the  BC  location  (Fig.  5a)  and  other  locations  along  the  axis  of  the  pancake, 
an  axisymmetric  texture  was  formed  for  the  primary  +  secondary  alpha  phases.  The 
textures  predicted  using  the  FEM  stresses  and  strain  increments  developed  at  810  °C  (Fig. 
6)  showed  that  specific  features  in  the  measured  pole  figures  could  be  explained  with  the 
variant  selection  rules  chosen  (Fig  5b  and  Fig  5c). 

In  the  case  for  which  all  variants  were  considered  possible  (Fig.  5b),  a  pole 
appeared  at  the  6  and  12  o’clock  positions  in  the  (1120)  pole  figure  and  a  fiber  along  the 
horizontal  equatorial  plane  of  the  (0002)  pole  figure.  Moreover,  fiber  textures  above  and 
below  the  horizontal  equatorial  planes  of  the  (0002)  and  (11 20)  pole  figures  were  also 
predicted.  These  features  were  also  observed  in  the  measured  pole  figures  (Fig.  5a). 

For  the  <11 1>  pencil-glide  variant-selection  approach,  which  allowed  several 
variants  to  form  from  a  single  beta  orientation  (Fig  5c)  the  fiber  textures  above  and  below 
the  horizontal  equatorial  planes  of  the  (0002)  and  (1120)  pole  figures  were  more 
pronounced.  Moreover,  two  strong  poles  appeared  diagonally  in  the  (0002)  and  (1120) 
pole  figures.  The  corresponding  measured  pole  figures  (Fig.  5a)  also  showed  poles  at 
these  locations,  suggesting  possibly  that  slip  system  activity  may  influence  which 
variants  will  form.  This  result  is  similar  to  that  found  in  the  work  of  Gey,  et  al.  [24],  who 
were  was  able  to  correlate  alpha-variant  selection  following  rolling  in  the  beta  phase  field 
to  the  most  active  {1 10}  <  1 1 1  >  and  {21 1}  <  1 1 1  >  slip  systems.  However,  the  lack  of 

predicted  poles  at  the  6  and  12  o’clock  positions  in  the  predicted  (11 20) pole  figure 
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indicates  that  this  variant-selection  approach  can  not  explain  a  number  of  the  features 
observed  in  the  measured  pole  figures  (Fig.  5a). 

A  (0002)  fiber-texture  component,  as  seen  by  the  poles  at  the  12  and  6  o’clock 
positions  in  the  (0002)  pole  figure  and  the  maximum  along  the  horizontal  equatorial 
plane  in  the  (1120)  pole  figure,  was  also  predicted  using  both  variant  selection  rules  (Fig 
5b  and  Fig  5c),.  Unfortunately,  these  types  of  texture  components  were  not  observed  in 
the  experimental  data  (Fig  5a).  Hence,  further  work  is  required  to  ascertain  the  origins  of 
this  texture  component. 

At  other  locations  in  the  pancake  forging  for  which  there  was  significant  rotation 
due  to  metal  flow,  such  as  RT,  a  more  complex  texture  was  formed  (Fig.  7a).  The 
stresses  and  strain  increments  developed  at  this  location  (Fig.  8)  were  dissimilar  to  those 
developed  at  the  BC  location  (Fig.  6).  Using  the  stresses  and  strains  developed  at  810  °C 
and  the  variant  selection  rules  based  on  <11 1>  pencil-glide  activity,  model  predictions 
(Fig  7c)  did  not  successfully  predict  the  majority  of  the  features  in  the  measured  pole 
figures  (Fig.  7a).  When  all  variants  were  considered  possible  (Fig  7b),  a  marked 
improvement  in  the  correlation  was  observed.  This  result  is  consistent  with  the 
microstructure  (Fig.  4d)  present  at  this  location,  in  that  somewhat  random  nucleation  and 
growth  of  the  alpha  phase  would  have  been  expected  to  have  occurred  in  the  production 
of  the  observed  basketweave  microstructure. 

The  cooling  rate  experienced  at  the  RT  location  was  sufficient  to  produce  the 
orthorhombic  (a”)  martensitic  phase  in  addition  to  the  secondary-alpha  platelets.  Based 
on  the  lattice  parameters  of  a”  (a  =  3.033,  b  =  4.924,  and  c  =  4.667  A)  deduced  in 
previous  x-ray  diffraction  work  [30,  31],  four  of  the  five  measured  pole  figures 
(((10  1  0),  (0002),  (10  1  1),  and  (1120))  for  the  primary  +  secondary  alpha  phase  were 
confounded  by  the  a”  martensitic  phase  due  to  diffraction-peak  overlap.  Moreover,  the 
a”  martensitic  phase  has  a 

(110)^  //  (001)a„ 

[in],//[ioiL„  (13) 

orientation  relationship  with  the  (3-phase  [32,  33],  with  perhaps  its  own  set  of  variant- 
selection  criteria  for  the  12  possible  variants.  Thus  the  a”  martentsitic  phase  will  be 
textured  in  some  manner  related  to  the  texture  of  the  p-phase  and  may  have  some  effect 
even  at  relatively  low  volume  fractions  on  the  calculated  orientation  distribution  function 
(ODF)  and  the  recalculated  pole  figures  for  the  primary  +  secondary  alpha  at  the  RT 
location  (Fig.  7a). 

From  a  quantitative  perspective,  the  model  predictions  of  the  alpha-phase  texture 
intensity  for  all  locations  were  reasonably  close  in  some  cases  and  high  in  others  for  the 
variant- selection  rules  examined.  At  the  BC  location  (i.e.,  Fig.  5),  the  predicted  texture 
using  all  variants  (Fig  5b)  was  comparable  to  the  measured  texture  (Fig  5a),  whereas  the 
modeled  texture  using  the  <11 1>  pencil-glide  activity  (Fig  5c)  was  4  to  5  times  that 
which  was  measured.  At  the  RT  location,  the  predicted  texture  strength  using  all  variants 
(Fig.  7b)  was  comparable  to  the  measurement  (Fig.  7a).  These  findings  can  be 
rationalized  in  part  due  to  the  fact  that  Taylor- type  models  typically  predict  textures  that 
are  stronger  than  those  measured.  Hence  in  regions  such  as  location  BC,  where  a  large 
amount  of  strain  is  accommodated  (Table  2),  an  over-prediction  in  the  texture  is  to  be 
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expected  and  will  be  compounded  by  the  use  of  any  type  of  variant-selection  biasing 
criterion  that  limits  which  variants  form  (Fig  5b).  In  addition,  the  number  of  discrete 
orientations  which  were  automatically  generated  to  represent  the  initial  billet  texture  of 
the  alpha  and  beta  phase  in  the  present  work  (-1000)  was  insufficient  to  fully  populate 
the  orientation  space.  As  a  result,  similar  to  a  coarse  grain  specimen  whose  texture  is 
measured  using  x-ray  diffraction,  the  simulated  pole  figures  were  spotty,  thus  limiting  the 
ability  to  make  quantitative  texture  predictions. 

Other  differences  in  the  measured-versus-predicted  pole  figures  of  a  quantitative 
nature  may  be  attributable  to  the  rapid  spatial  variation  in  the  cooling  profiles  in 
specimens  close  to  the  outer  extremities  of  the  forging.  In  specimens  close  to  the 
pancake  outer  surface,  such  as  the  RT  specimen,  the  modeled  stresses,  cooling  rates,  and 
microstructures  varied  significantly  over  the  surface  of  the  specimen  in  spite  of  the  large 
initial  billet  and  final  pancake  forging  dimensions.  Moreover,  the  modeled  cooling  rates 
and  observed  changes  in  the  microstructures  developed  were  consistent  with  the 
formation  of  a  low  volume  fraction  of  the  textured  a”  martensitic  phase,  which  would 
also  contribute  to  the  observed  quantitative  differences  in  the  pole  figures. 

For  specimens  at  the  interior  of  the  pancake  forging,  such  as  the  BC  location, 
cooling-rate  gradients  were  present,  but  to  a  much  smaller  degree.  As  a  result,  the 
microstructures  were  relatively  homogeneous  over  the  entire  sample  surface,  and  the 
gradients  were  concluded  to  play  a  minor  role  in  the  differences  in  the  intensities  of  the 
various  components  of  the  measured  and  predicted  pole  figures. 

A  comparison  of  the  measured  textures  for  RT  samples  that  were  solution  treated 
and  furnace  cooled  (Fig.  9a)  (to  retain  solely  the  equiaxed  alpha  phase  and  its  associated 
deformation  texture)  or  water  quenched  (Fig.  7a)  showed  significant  differences,  thus 
illustrating  the  contribution  of  the  transformation  texture  to  the  overall  alpha-phase 
texture.  Furthermore,  the  deformation-texture  prediction  for  the  furnace-cooled  RT 
sample  showed  qualitative  agreement  with  the  measurement  (Figure  9a  versus  Figure  9b), 
but  the  texture  intensity  was  again  approximately  twice  as  strong.  Because  a  Taylor-type 
crystal  plasticity  model  was  used,  such  a  result  was  again  not  unexpected. 

6.  Summary  and  Conclusions 

The  development  of  crystallographic  texture  during  the  forging  and  heat  treatment 
of  two-phase  (alpha/beta)  titanium  alloys  has  been  modeled  by  coupling  software  to  treat 
strain  partitioning  between  the  phases,  crystal  rotations  due  to  metal  flow  and  slip,  and 
variant  selection  during  the  decomposition  of  beta  to  form  secondary  (platelet)  alpha. 
The  method  was  evaluated  using  an  industrial-scale  Ti-6A1-4V  pancake  forging.  A 
comparison  of  measured  and  predicted  textures  showed  qualitative  and  in  some  cases 
quantitative  agreement.  However,  several  refinements  are  needed  such  as  improvement 
in  the  deformation  texture  model  to  reduce  the  tendency  for  over-prediction  of  texture 
intensity  and  methods  to  account  for  possibly  competitive  variant-selection  mechanisms. 
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Table  1.  Subset  ofcubic  symmetry  operators  used  to  transform  the  beta  phase  into 
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Table  2.  Dependence  of  the  volume  fraction  of  primary  alpha 
and  effective  strain  on  location  in  the  pancake  forging 


Location 

Volume  Fraction  of  Primary  Alpha  (%) 

Effective  Strain 

BC 

46 

2.635 

MC 

48 

2.011 

RC 

42 

1.266 

BT 

38 

0.055 

MT 

39 

0.599 

RT 

30 

1.105 
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Figure  Captions 

Fig.  1.  Backscattered-electron  image  of  microstructure  developed  in  Ti-6A1-4V  at  955°C 
followed  by  water  quenching.  The  dark,  equiaxed  phase  is  primary  alpha.  The 
gray,  transformed  phase  was  beta  at  the  heat-treatment  temperature. 

Fig.  2.  Specimens  extracted  from  the  Ti-6A1-4V  pancake  forging  for  microstructure  and 
texture  examination.  BC:  bore-midplane,  MC:  midradius-midplane,  RC:  rim- 
midplane,  BT:  bore-top,  MT:  midradius-top,  and  RT:  rim-top  locations. 

Fig.  3.  FEM-model  predictions  of  temperature  transients  at  different  locations  in  the 
pancake  forging  during  quenching  following  final  heat  treatment. 

Fig.  4.  Backscattered-electron  images  of  microstructures  developed  in  the  Ti-6A1-4V 
pancake  forging  following  water  quenching  from  the  final  heat  treatment  at  the  (a) 
bore-midplane,  (b)  rim-midplane,  (c)  bore -top,  and  (d)  rim-top  locations. 

Fig.  5.  Alpha-phase  pole  figures  for  the  bore-midplane  (BC)  location  following  forging 
and  heat  treatment:  (a)  measured,  (b)  simulated  with  all  variants,  (c)  simulated 
with  preferential  variant  selection  based  on  <1 1 1>  pencil-glide. 

Fig.  6.  FEM  predictions  at  the  bore-midplane  (BC)  location:  (a)  stress  and  (b)  strain 
components  and  temperature  as  a  function  of  time  during  water  quenching 
following  solution  heat  treatment  of  the  forging. 

Fig.  7.  Alpha-phase  pole  figures  for  the  rim-top  (RT)  location  following  forging  and  heat 
treatment:  (a)  measured,  (b)  simulated  with  all  variants,  (c)  simulated  with 
preferential  variant  selection  based  on  <1 1 1>  pencil-glide. 

Fig.  8.  FEM  predictions  at  the  rim-top  (RT)  location:  (a)  stress  and  (b)  strain  components 
and  temperature  as  a  function  of  time  during  water  quenching  following  solution 
heat  treatment  of  the  forging. 

Fig.  9.  Primary-alpha  pole  figures  for  the  rim-top  (RT)  location  following  forging  + 
solution  heat  treatment  +  furnace  cooling:  (a)  measured  and  (b)  simulated. 
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